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ABSTRACT 

We present a measurement of the spatial clustering of submillimetre galaxies (SMGs) at 
z = 1-3. Using data from the 870 /im LABOCA submillimetre survey of the Extended Chan- 
dra Deep Field South, we employ a novel technique to measure the cross-correlation between 
SMGs and galaxies, accounting for the full probability distributions for photometric redshifts 
^y-^ I of the galaxies. From the observed projected two-point cross-correlation function we derive 

■ the linear bias and characteristic dark matter halo masses for the SMGs. We detect cluster- 

ing in the cross-correlation between SMGs and galaxies at the > 4(7 level. Accounting for 
the clustering of galaxies from their autocorrelation function, we estimate an autocorrelation 
length for SMGs of tq = 7.7^2 3 ^P'^ assuming a power-law slope j — 1.8, and derive 
a corresponding dark matter halo mass of log(Mhaio[^"^ Mq]) = 12.8t°j. Based on the 
evolution of dark matter haloes derived from simulations, we show that that the 2 = descen- 
dants of SMGs are typically massive (~2-3 L*) elliptical galaxies residing in moderate- to 
, high-mass groups (log(Afhaio[/i~"'^ -^©l) = 13. Sj^Q'g). From the observed clustering we esti- 

^ • mate an SMG lifetime of 100 Myr, consistent with lifetimes derived from gas consumption 

\ times and star-formation timescales, although with considerable uncertainties. The clustering 

of SMGs at z ^ 2 is consistent with measurements for optically-selected quasi-stellar ob- 
jects (QSOs), supporting evolutionary scenarios in which powerful starbursts and QSOs oc- 
cur in the same systems. Given that SMGs reside in haloes of characteristic mass 6 x 10^^ 
Mq, we demonstrate that the redshift distribution of SMGs can be described remarkably 
well by the combination of two effects: the cosmological growth of structure and the evolution 
of the molecular gas fraction in galaxies. We conclude that the powerful starbursts in SMGs 
likely represent a short-lived but universal phase in massive galaxy evolution, associated with 
the transition between cold gas-rich, star-forming galaxies and passively evolving systems. 

Key words: galaxies: evolution - galaxies: high-redshift - galaxies: starburst - large-scale 
structure of the Universe - submillimetre. 
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1 INTRODUCTION 

Submillimetre galaxies (SMGs) are a population of high-redshift 
ultraluminous infrared galaxies (ULlRGs) selected through their 
reds hifted far-infrared emission in the submillimetre waveband 
(e.g. ISmail. Ivisori & Blainll997l;lBM-ger et alJ[T998l ; lHughes et"ai] 
Il998l ; lBlainetalJl2002h . The redshift dist ribution of this popu- 
lation appears to peak at z ~ 2.5 (e.g., IChapmanet"ani2003l . 
I2OO5I : IWardlow et al.ll2oTTh . so that SMGs are at their common- 
est around the same epoch as the peak in powerful active galactic 
nuclei (AGN) and specifically quasi-s tellar objects (QSOs) (e.g., 
iRichards et ZII2OO6I : IXssef et al.ll201 ih . This correspondence may 
indicate an evolutionary link between SMGs and QSOs, similar 
to th at suggested at low redshift between ULIRGs and QSOs by 
I Sanders et ah tl988l . However there is little direct overlap (~a 



few percent) between the hieh-redshift SMG and QSO popula- 
tions (e.g.. Page et al. 2004; Chapman et al. 2005: IStevens et alj 
[2005; Alexander et al. 2008; Wardlow et al. 2011). The immense 
far-infrared luminosities of SMGs are widely believed t o arise from 
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Alexander et aill2005l : IPope et alj|2008l ; iTacconi et al.ll2006l 
Ivison et al.ll201ir) . suggesting that they may represent the 



formation phase of the most massive local galax ies: giant ellipti- 
cals (e.g.. lEales et al.lll999l ; [Swinbank et alhood) . 

SMGs and QSOs may thus represent phases in an evolu- 
tionary sequence that eventually results in the population of lo- 
cal massive elliptical galaxies. This is a compelling picture, but 
testing the evolutionary links is challenging due to the lack of an 
easily-measured and conserved observable to tie the various pop- 
ulations together. For example, the stellar masses of both QSOs 
and SMGs are difficult to measure reliably due to either the 
brigh t ness of the nuclear em ission in the QSOs (e.g.. lCroom et alj 
|2004 iKotilainen et al.ll200^ or strong dust obscuration and po- 
tentially comp lex star-formation histories for the SM Gs (e.g., 
Hainline et al. 201 1 ; 'Wardlo w et al.l [20I1I : but see also iDunlopI 



2011 ; Michalowski et al„ ^QJUD. while the details of the high- 
redshift star formation that produced loc al massive elliptical g alax- 
ies are likewise poorly constrained (e.g.. lAllanson et alj2009l) . De- 
riving dynamical masses for QSO hosts from rest-frame optical 
spectroscopy is difficult due to the very broad emission lines from 
the AGN, while dynamical mass measurements using CO emis- 
sion in gas-rich QSOs are also challenging, due to the potential 
non-isotropic orientation of the QSO hosts on the sky and the 
lack of high-resolu tion velocity fields necessary to solve for this 
dCoDpin et al.H2008h . as well as the general difficulties in model- 
ing CO kinematic s (e.g.. lTacconi et alJl2006l ; lBothwell et alj|20ld : 
lEngeletalJl2010l) . 

Another possibility is to compare source populations via 
the masses of their central black holes. For QSOs and the pop- 
ulation of SMGs that contain broad-line AGN, the black hole 
mass can be estimated using virial techniques based on the 
broa d emission lines (e. g., Vestergaard 2002; Peterson et al. 20041; 
Vestergaard & Peterson! Bood ; iKollmeier et alj |2006| ; IShen et alj 
20081) . Such studies generally find that SMGs have small black 



holes relative to the local black hole-gal axy mass relations (e.g., 
[Alexander et alj2008l ; ICMTera et alj201 ih . while the black holes in 
2 ~ 2 QSOs tend to lie above the local relatio n, with masses sim- 
ilar to those in local massive elliptical s (e.g., iDecarli et al]|2010l ; 
iBennert et alj2010l ; lMerloni et all2O10l) . These results suggest that 
SMGs represent an earlier evolutionary stage, prior to the QSO 
phase in which the black hole reaches its final mass. However, high- 
redshift virial black hole mass estimates are highly uncertain (e.g.. 
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Figure 1. Two-dimensional distribution of the 50 LESS SMGs and 
^ 50,000 IRAC galaxies in the ECDFS that are used in our analysis. The 
SMG s shown repres ent the subset of the 126 SMGs in the full LESS sam- 
ple dWeifi et al.ll2009n that are in the redshift range 1 < 2 < 3 and are 
in regions of good photometry, and so aie used in this analysis. The IRAC 
galaxies are chosen to reside at 0.5 < z < 3.5. The SMGs are shown here 
individually, while the density of galaxies is given by the grayscale. The 
blank areas represent regions which are excluded from the analysis, includ- 
ing areas of poor photometry (for example around bright stars) or additional 
sources identitied by eye in the vicinity of SMG, as discussed in [|2] The 
high density of IRAC galaxies in the field enables an accurate measurement 
of the SMG-galaxy cross-correlation function. 
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Figure 2. Redshift distributions for the IRAC galaxy sample in the redshift 
range 0.5 < 2 < 3.5 (dotted line), and the SMG sample in the range 
1 < 2 < 3 (solid line). The histogram for galaxies has been scaled so 
that the distribution can be directly compared to that of the SMGs. Also 
shown is the redshift distribution for 11,241 galaxies (dashed line) selected 
to match the overlap in the redshift distributions of the SMGs and galaxies, 
as used in the galaxy autocorrelation measurement ( ^3.2t . For the SMGs, 
44% have spectroscopic redshifts, while the remainder of the SMGs and all 
the IRAC galaxies have redshift estimates from photometric redshift calcu- 
lations ( Wardlow et al. 2011.) . 
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iMarconi et al] l2008l : iFine et alj l20ld iNetzer & Marzianil boiot) 
and may suffer from significant selection effects (e.g., Lauer et all 
[2007; Siien& Kelly 2010; Kelly et al. 2010), and so conclusions 
about connections between populations are necessarily limited. 

The difficulties discussed above lead us to take another route 
to compare SMGs to high-redshift QSOs and low-redshift el- 
lipticals: through their clustering. Spatial correlation measure- 
ments provide information about the characteristic bias and hence 
mass of the hgJoes in which galaxies reside (e.g.. Kaiser 1984; 
iBardeen et al ]|1986'), and so provide a robust mass estimate that is 
free of many of the systematics in measuring stellar or black hole 
masses. The observed clustering of SMGs and QSOs can thus allow 
us to test whether these populations are found in similar haloes and 
so may evolve into each other over short timescale s. With knowl- 
edge of how haloes evolve over cosmic time (e.g.. lLaceY & Colj 
ll993l ; lFakhouri, Ma & Bovlan-Kolchinll201ol). we can also exp lore 



the links to modem eUiptical galaxies (e.g., Overzier et alj2003l) . as 
well as the higher-redshift progenitors of SMGs. Clustering mea- 
surements can also provide constraints on theoretical studies that 
explore the nature of SMGs in a cosmological context. Recent mod- 
els for SMGs as relatively long-lived (> 0.5 Gyr) star formation 
episodes in the most ma ssive galaxies, d riven by the early collapse 
of the dark matter halo iXia et al. 201 1), or powered by steady ac- 
cretion of intergalactic gas (Dave et al. 2010), yield strong cluster- 
ing for bright sources (850 pm fluxes > a few mjy) with correlation 
lengths ro > 10 Mpc. In contrast, models in which SMGs are 
short-lived bursts in less massive galaxies, with large luminosities 
produced by a t 
weaker clustering 

[mil). 

Attempts to measure the clustering of SMGs from their 
projected two-dimensional distri bution on the s k y have for 
the most part been ambiguous |Scott et alj |2002|; [Borys et alj 



I top-heavy initial mass fun ction, predic t significa ntly 
ring with ro ~ 6h~^ Mpc iAlmeida, Baugh & LacevI 



2003'; 'Webb et al.' ^2003 



WeiB et al.' '20091; Iwilliams et alj 1201 ll ; 

Lindner et al. 2011). WeiB et al. (2009) used the largest, contigu- 
ous extragalactic 870-/im survey (of the Extended Chandra Deep 
Field South; ECDFS), to derive the clustering of ^ 5-mJy SMGs 
from their projected distribution on the sky. They e stimated a cor- 
relatio n length of 13 ± 6h^^ Mpc. Most recently. Iwilliams et alj 
( I2OIII) analysed a 1100-^im survey of a region of the COSMOS 
field and placed l-cr upper limits on the clustering of bright SMGs 
(with apparent 870-/im fluxes > 8-10 mJy) of > 6-12 Mpc. 

Other work has attempted to improve on angular correla- 
tion measurements by including redshift information. Using the 
spectroscopic redshift survey of 73 SMGs with 870-Atm fluxes o f 
5 mJy spread a cross seven fields from Ichapman et alj j2005l) . 
iBlain et ^ ( l2004h estimated a clustering amplitude from the num- 
bers of pairs of SMGs within a 1000-kms^^ wide velocity win- 
dow. They derived an effective correlation length of 6.9 ± 2.1 
/i~^Mpc, suggesting that SMGs are strongly clust ered. How- 
ever t heir methodology was subsequently criticised by jAdelbergej 
l2005l) . who suggested that accounting for angular clustering of 
sources and the redshift selection function significantly increases 
the uncertainties. Using data from the Chandra Deep Field-North, 
lBlakeetalJl2003) computed the angular cross-correlation between 
SMGs and galaxies in slices of spectroscopic and photometric red- 
shift. They obtained a significant SMG-galaxy cross-correlation 
signal, with hints that SMGs are more strongly clustered than the 
optically-selected galaxies, although with only marginal (~2a) 
significance. Previous work has therefore pointed toward SMGs 
being a strongly clustered population, but their precise clustering 



amplitude, along with their relationship to QSOs and ellipticals, re- 
mains uncertain. 

To make improved measurements of th e clustering of SM Gs, 
we need either much larger survey areas (see lCoorav et al.l201(ll for 
a wide-field clustering measurement for far-IR detected sources) or 
the inclusion of redshift information (to allow us to reduce the ef- 
fects of projection o n our cluste ring measurements). To this end, we 
have reanalysed the lWeifi et alj l2009i) survey of ECDFS using new 
spectroscopic a nd photometric redsh ift constraints on the counter- 
parts to SMGs jWardlow et alj[2oTll) as well as a large catalogue 
of "normal" (less-active) galaxies in t he same field. We employ a 
new clustering analysis methodology l lMvers. White & Ball 2009h 
to calculate the projected spatial cross-correlation between SMGs 
and galaxies, to obtain the tightest constraint to date on the cluster- 
ing amplitude of SMGs. 

This paper is organised as follows. In § 2 we introduce the 
SMG and galaxy samples, and in § 3 we give an overview of the 
methodology used to measure correlation functions and estimate 
dark matter (DM) halo masses. In § 4 we present the results, explore 
the effects of photometric redshift errors, compare with previous 
measurements, and discuss our results in the context of the physical 
drivers, lifetimes, and evolutionary paths of SMGs. In § 5 we sum- 
marise our conclusions. Throughout this paper we assume a cos- 
mology with Qui = 0.3 and Qa = 0.7. For direct comparison with 
other works, we assume Ho = 70 km s~^ Mpc~^ (except for co- 
moving distances and DM halo masses, which are explicitly given 
in terms of h — Ho / ( 1 00 km s ~ ^ Mpc ~ ^ ) ) . In order to easily com- 
pare to estimated halo m asses i n other recent works on QSO clus- 
tering (e.g., Croom e t alj|2005l : iMvers et aLll2006l ; Ida Angela et alj 
12003 ; iRossetal J 1200^ we assume a normalisation for the mat- 
ter power spectrum of as = 0.84. All quoted uncertainties are la 
(68% confidence). 



2 SMG AND GALAXY SAMPLES 

Our SMG sample comes from the survey of the ECDFS using the 
Large APEX BOlometer CAmera dSiringo et alll2009. LA B OCA) 
on the Atacama Pathfinder Experiment dOtisten et alj200d APEX) 
12-m telescope (the LA BOCA ECDFS Submillimetre Survey, or 
LESS; IWeifi etai]|2009l) . LESS mapped the full 0.35 deg^ ECDFS 
to a 870-^m noise level of 1 .2 mJy beam~^ a nd detected 126 
SMGs at > 3.7a significance dWeiB et al]|2009i equivalent to a 
false-detection rate of ~ 4%). Radio a nd mid-infr a red co unter- 
parts to LESS SMGs were identified by 'Biggs et all ( 1201 ll) using 
a maximum-likelihood technique. Spectroscopic and photometric 
redshift s were obtained for a significant fraction of these counter- 
parts bv lWardlowetalJ ( [20T l|) and we refer the reader to that work 
for more details. For this study, we restrict our analysis to the 50 
SMGs that have secure counterparts at z = 1-3 and do not lie close 
to bright stars (as discussed below). The upper limit of 2: = 3 on the 
sample is included to maximize overlap in redshift space with the 
galaxy sample, in order to obtain a significant cross-correlation sig- 
nal, while the lower bound of 2; = 1 is included to prevent the SMG 
sample from being biased toward low redshifts. Of the SMGs in the 
sample, 22 SMGs (44%) have spectroscopic redshifts (Danielson 
et al., in preparation) and the remainder have photometric redshift s 
with a typical precision of CTz/(1 + 2) ~ 0.1 jWardlow et alj2oTl[) . 
The 8 70-/jm flux distrib ution for the SMGs having secure counter- 
parts (Biggs et al .ll20Ilh is consistent with that for all LESS SMGs 
I WeiB et aL (.200% . indicating that the requirement that SMGs have 
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secure counterparts does not strongly bias the fluxes of our SMG 
sample. 

For the cross-correlation analysis, we also require a compar- 
ison population in the same field. For this we adopt the ~ 50,000 
galaxies detected in the Spitzer IR AC/MUS YC Public Legacy Sur- 
vey in the Extended CDF- South jPamen et al.ll201lh . We use an 
IRAC selected sample to ensure that each galaxy has photom- 
etry in a sufficient number of bands, and over a wide enough 
wavelength range, to allow robust estimates of photometric red- 
shift. Photo-zs are calculated using template fits to the optical and 
IRA C photo metry in an identical method to that used for the SMGs 
(see Ward low et al. 2011). The fits are performed with HYPER-Z 
jBolzonena , Miralles & Pellg 200Q) and the resulting redshift dis- 
tribution, compared to that for the SMGs, is shown in Figure [2] 
The photometric analysis uses chi-squared minimisation, which al- 
lows the calculation of confidence intervals for the best-fit redshift. 
These can be presented as a probability distribution function (PDF) 
for the redshift, or equivalently, the comoving line-of-sight distance 
X (calculated for our assumed cosmology). We define the PDF for 
each galaxy as /(x), where J f{x)dx ~ 1- Examples of the PDFs 
for the galaxies are shown in Figure [3] 

Finally, in order to calculate the correlation functions, we 
first create random catalogues of "galaxies" at random positions 
within the actual spatial coverage of our survey. Like many fields, 
the ECDFS contains several bright stars with large haloes, around 
which few galaxies are detected. Th erefore, we use the backg round 
map produced by S EXTRACTOR jSertin & AmoutsI \l99Q) from 
the combined IRAC image during the source extraction proce- 
dure to create a mask. This mask is applied to the random cat- 
alogues, the SMGs and the IRAC galaxies, so that the positions 
of the random galaxies are unbiased with respect to the SMG 
and IRAC galaxy samples, and thus the mask does not affect the 
cross-correlation measurement. As discussed in iBiggs et aljfeoilh 
and I Wardlow etai] §011), some of the SMG identifications were 
performed manually by examining the regions around the SMGs. 
These additional sources are excluded from the clustering analysis 
so as not to bias the results. The sky positions of the SMGs and 
galaxies that are outside the masked regions are shown in Figure[T] 



3 CORRELATION ANALYSIS 

To measure the spatial clustering of SMGs, we can in principle 
derive the autocorrelation of the SMGs themselves. However, as 
we have discussed, current SMG samples are too limited in size and 
available redshift information to make this feasible. Alternatively, 
we can measure the croM-correlation of a population with a sample 
of other sources (for ex ample, less-active ga l axies) which populate 
the same volume (e.^.. iGawiser et ^ 2001 : Adelberger & Steidel 



l2005l : lBlake et alj|2006| - 1coil et al.ll2007l : iHickox et alj|2009l) . The 

much larger number of galaxies in the ECDFS (~ 1000 x more 
than the SMGs in a comparable redshift range) allows far greater 
statistical accuracy in the measurement of clustering. 

To calculate the real-space projected cross-correlation func- 
tion Wj,( R) between SMGs and galax ies we employ a method de- 
rived bv'Mvers. White & Balll ( l2009 l). This method enables us to 
take advantage of the full photo-z PDF for each galaxy, by weight- 
ing pairs of SMGs and galaxies based on the probability of their 
overlap in redshift space. This method allows us to calculate the 
SMG-galaxy cross-correlation using the full sample of 2 ~ 50, 000 
IRAC galaxies, while the derive the clustering of the galaxies them- 
selves using a smaller sample that is selected to match the overlap 
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Figure 3. Example probability distribution functions for thi'ee IRAC galax- 
ies and an SMG. We mark the "best" (peak) comoving distance for each 
galaxy. Note that for each galaxy in this example, the line-of-sight distance 
between the "peak" redshift of the galaxy and the SMG redshift is far too 
large for them to be physically associated. However, because of the uncer- 
tainty in the galaxy redshifts (shown by the PDFs), there is a non-negligible 
probability that the galaxies lie close to the line-of-sight distance of the 
SMG. 



in the redshift distributions of the galaxies and SMGs. Our cluster- 
ing analysis is identical in mo st respects to t he QS O-galaxy cross- 
correlation study presented in lHickox et alj ^201 ll hereafter HI 1). 
Because the method is somewhat i nvolv ed, we present only the key 
details here and refer the reader tolHllI for a full discussion. 



3.1 Cross-correlation method 

The two-point correlation function ^(r) is defined as the probability 
above Poisson of finding a galaxy in a volume element dV at a 
physical separation r from another randomly chosen galaxy, such 
that 



dP = n[l+^{r)]dV, 



(1) 



where n is the mean space density of the galaxies in the sample. 
The projected correlation function 'Wp{R) is defined as the integral 
of 5(r) along the line of sight. 



Wp{R) 



^{R, TT)d-K, 



(2) 



where R and tt are the projected comoving separations between 
galaxies in the directions perpendicular and parallel, respectively, 
to the mean line of sight from the observer to the two galaxies. 
By integrating along the line of sight, we eliminate redshift-space 
distortions owing to the peculiar motions of galaxies, which dis- 
tort the line-of-sight distances measured from redshifts. Wp{R) has 
been used to measure correlations in a number of surveys (e.g. , 
'Zehavi et al."2005'; Li et al. 2006; Gilli et al. 2007; Coil et al. 200 
2008; Wake et al 2008a; Myers, White & Ball 2009.; Hickox et al. 



20095 ICoil et al l lToogI; iGilli et alj |2009|; iKrumpe, Mivaii & Coi 



2010; Donoso et alj2010l ; lHickox et alj201ll ; IStarikova et alj201ll ; 
AUevato et al. 201j|). 

In the range of separations 0.3 ^ r < 50 Mpc, ,^(r) for 
galaxies and QSOs is roughly observed to be a power-law. 



C(r) = (r/ro)" 



(3) 
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with 7 typically 1.8 (e.g., IZehavi et alj l2005l : ICoil etalj[20o3. 
1 20071 : IRoss et al]|200 9V For sufficiently large TTmax such that we 
average over all line-of-sight peculiar velocities, 'Wp{R) can be di- 
rectly related to ^(r) (for a power law parameterisation) by 



r(i/2)r[(7 - 1)/2] 
r(7/2) 



(4) 



To calculate Wp{R) for the cro ss-cor relation between SMGs 
and galaxies, we use the method of lMOSL which accounts for the 
photometric redshift probability distribution for each galaxy indi- 
vidually. Following M09, the projected cross-correlation function 
can be calculated using: 



Wp(R) 



NrNs ^ C: 



DsDajR) 
DsRg{R) 



id 



where 

'^'d ~ fid / ^ ^ fid 



(5) 



(6) 



Here R is the projected comoving distance from each SMG, for 
a given angular separation 6 and radial comoving distance to the 
SMG of X*, such that R = x*^- DsDg and DsRg are the 
number of SMG-galaxy and SMG-random pairs in each bin of 
R, and Ns and Nr are the total numbers of SMGs and random 
galaxies, respectively. is defined as the average value of the ra- 
dial PDF f{x) for each galaxy i, in a window of size Ax around 
the comoving distance to each spectroscopic source j. We use 
Ax = 100 /i~^Mpc to effectively eliminate redshift space dis- 
tortions, although the results are insensitive to the details of this 
choice. We refer the reader to M09 and HI 1 for a detailed deriva- 
tion and discussion of these equations. In this calculation as well 
as in the galaxy autocorrelation, we account for the integral con- 
straint as described in lHlll . This correction increases the observed 
clustering amplitude by «15%. 



3.2 Galaxy autocorrelation 

To estimate DM halo masses for the SMGs, we calculate the rel- 
ative bias between SMGs and galaxies, from which we derive the 
absolute bias of the SMGs relative to DM. As discussed below, 
calculation of absolute bias (and thus halo mass) requires a mea- 
surement of the autocorrelation function of the IRAC galaxies. The 
large size of the galaxy sample enables us to derive the clustering 
of the galaxies accurately from the angular autocorrelation function 
io{9) alone. Although we expect the photometric redshifts for the 
IRAC galaxies to be reasonably well-constrained (as discussed in 
§[2j, by using the angular correlation function we minimize any un- 
certainties relating to individual galaxy photo-zs for this part of the 
analysis. The resulting clustering measured for the galaxies has sig- 
nificantly smaller uncertainties than that for the SMG-galaxy cross- 
correlation. 

We calculate the angu lar autocorrelation function ijj{9) using 
the lLandv & Szalavl ( Il993h estimator: 

uj{e)^-^{DD-2DR + RR), (7) 

where DD, DR, and RR are the number of data-data, data- 
random, and random-random galaxy pairs, respectively, at a sep- 
aration 6, where each term is scaled according to the total numbers 
of SMGs, galaxies, and randoms. 

The galaxy autocorrelation varies with redshift, owing to the 



evolution of large scale structure, and because the use of a flux- 
limited sample means we select more luminous galaxies at higher 
z. This will affect the measurements of relative bias between SMGs 
and galaxies, since the redshift distribution of the SMGs peaks at 
higher z than that for the galaxies and so relatively higher-z galax- 
ies dominate the cross-correlation signal. To account for this in our 
measurement of galaxy autocorrelation, we randomly select galax- 
ies based on the overlap of the PDFs with the SMGs in comoving 
distance (in the formalism of §[3T|this is /ij for each galaxy, av- 
eraged all SMGs). We select the galaxies so their distribution in 
redshift is equivalent to the weighted distribution for all galaxies 
(weighted by {fij})- The redshift distribution of this galaxy sam- 
ple is shown in Figure |2l We use this smaller galaxy sample to 
calculate the angular autocorrelation of IRAC galaxies. 



3.3 Uncertainties and model fits 

We estimate uncertainties on the clustering directly from the data 
using bootstrap resampling. Following HI 1, we divide the field into 
a small number of sub-areas (we choose A'sub = 8), and for each 
bootstrap sample we randomly draw a total of SN^uh sub-areas 
(with replacement), which has been shown to best a pproximate the 
intrinsic uncertainties in the clustering amplitude jNorberg et alj 
2009). To account for shot noise owing to the relatively small size 
of the SMG sample, we take the sets of SN^uh bootstrap sub- 
areas and randomly draw from them (with replacement) a sample 
of sources (SMGs or galaxies) equal in size to the parent sample; 
only pairs including these sources are used in the resulting cross- 
correlation calculation. We use the bootstrap results to derive the 
covariance between differen t bins of R, calculating the covariance 
matrix using Equation 12 of lHl ll . 

We fit the observed Wp(R) with two models: a power law and 
a simple bias model (described in ij |3.4| (. We compute model param- 
eters by minimising (taking into account the covariance matrix 
as in Equation 13 of iHUh and derive la errors in each parame- 
ter by the range for which Ax'^ = 1. We use the same formal- 
ism for computing fits to the angular correlation functions, where 
<-o{6) = A9^^ . We convert A and <5 to real-space clustering param- 
eters ro and 7 following the procedure described in § 4.6 of iHl U 



3.4 Absolute bias and dark matter halo mass 

The masses of the DM haloes in which galaxies and SMGs reside 
are reflected in their absolute clustering bias fcabs relative to the DM 
distribution. The linear bias 6^i,s is given by the ratio of the autocor- 
relation function of the galaxies (or SMGs) to that of th e DM . We 
determine feabs following the method outlined in § 4.7 of lHl iL sim- 
ilar to the approach used previously by a number of studies (e.g., 
iMve rs et al.l200^,l2007l : ICoil et alj2007ll2008l , l2009l : rHickox et alj 
2009); in what follows we briefly describe this procedure. 

We first calculate the two-point autocorrelation of DM as 
a fun ction of redshift. We use the HALOFIT code of ISmith et alj 
( |2003[) assuming our standard cosmology, and the slope of the ini- 
tial fluctuation power spectrum, F — Q.mh = 0.21, to derive 
the DM power spectrum, and thus its projected correlation func- 
tion 10°^ (/?), averaged over the redshift distribution for which the 
SMGs and galaxies overlap. We then fit the observed Wp (R) of the 
SMG-galaxy cross-correlation, on scales 0.3-15 Mpc, with a 
model comprising a simple linear scaling of w^^{R). The best-fit 
linear scaling of the DM correlation function corresponds to bsbc, 
the product of the linear biases for the SMGs and galaxies, respec- 
tively. This simple model produces a goodness-of-fit comparable 
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Figure 4. The projected SMG-galaxy cross-coiTelation function (derived 
using Equation|5j. Uncertainties are estimated from bootstrap resampling. 
A power-law fit to Wp{R) is shown by the solid line, and the projected 
correlation function for DM is shown by the dotted line. Fits are performed 
over the range in separation of R = 0.3-15 h^^ Mpc. Both the power law 
model with 7 = 1.8 and a linear scaling of the DM correlation function 
provide satisfactory fits to the observed Wp{R). Together with the observed 
galaxy autocorrelation, this measurement yields the clustering amplitude 
and DM halo mass for the SMGs, as desciibed in §|4] 



Figure 5. The angular autocorrelation function of IRAC galaxies, selected 
to match the overlap of the SMGs and galaxies in redshift space. Uncer- 
tainties are estimated from bootstrap resampling. The angular correlation 
function for DM, evaluated for the redshift distributions of the galaxies, is 
shown by the dotted gray line. The power law fit was performed on scales 
0.3'-10' and is shown as the solid line. Both the power law model with 
5 = 0.8 and a linear scaling of the DM correlation function provide sat- 
isfactory fits to the observed uj{6). The observed amplitude of the galaxy 
autocon'elation yields the absolute bias of the galaxies, which we use to 
obtain the absolute bias and DM halo mass of the SMGs. 



to that of the power-law model in which the slope 7 is allowed to 
float. 

To determine bs we therefore need to estimate ba- We obtain 
be for the galaxies from their angular autocorrelation in a similar 
manner to that applied to the SMG-galaxy cross-correlation. Again 
we calculate the autocorrelation for the DM uidm{0), by integrat- 
ing the power spe ctrum from HALOFIT using Equation (A6) of 
iMvers et al.l jZOOTh . We fit the observed u){9) with a linear scaling 
of uiDM{d) on scales 0.3'-10' (corresponding to 0.3-10 Mpc 
at 2; = 2). This linear scaling corresponds to feg and thus (combined 
with the cross-correlation measurement) yields the SMG bias bs- 
Finally, we convert bp and bg to Mhaio usi ng the prescription of 
ISheth. Mo & Tormed JzOOlh . as described in lHlll . This character- 
istic Mhaio corresponds to the top-hat virial mass (see e.g.. |Peeblej 
I1993I and references therein), in the simplified case in which all 
objects in a given sample reside in haloes of the same mass. This 
assumption is justified by the fact (as discussed below in ii l4.4t that 
SMGs have a very small number density compared to the popu- 
lation of similarly-clustered DM haloes, such that it is reasonable 
that SMGs may occupy haloes in a relatively narrow range in mass. 
We note that this method differs from some prescriptions in the lit- 
erature which assume that sources occupy all haloes above some 
minimum mass; this is particularly relevant for populations with 
high number densities that could exceed the numbers of available 
DM haloes over a limited mass rang e. Given the halo mass func- 
tion at 2 ~ 2 (e.g., iTinlcer et alj[2008i) the derived minimum mass 
is typically a factor of ~2 lower, for the same clustering amplitude, 
than the "average" mass quoted here. 



4 RESULTS AND DISCUSSION 

The projected cross-correlation function of the SMG sample with 
the IRAC galaxies is shown in Figure|4] We plot the best-fit power- 



law model, and show the correlation function of the DM calcu- 
lated as in § 13.41 which we fit to the data through a linear scal- 
ing. The power-law and linear bias fit parameters are presented 
in in Table [T] For SMGs the observed real-space projected cross- 
correlation is well-detected on all scales from 0.1-15 Mpc, and 
the power-law fits return 7 1.8, similar to m any previous corre- 
lation function m easurements fo r galaxies (e .g.,'Zehavi et al.*20()j; 
ICoil et alj200il and QSOs (e.g.. lCoil et Zr2 007: Ross et al. 200^. 
The best-fit parameters for the SMG-galaxy cross-correlation are 
ro,SG = 5.3 ± 0.8 h'^ Mpc, 7 = 1.7 ± 0.2. If we fix the value 
of 7 to 1.8, we obtain ro,sG ~ 5.1 ± 0.6 Mpc, corresponding 
to a clustering signal that is significant at the > 4a level, the most 
significant measurement of SMG clustering to date. From the fit of 
the DM model, we obtain bsba = 5.83 ± 1.36. 

We next compute the autocorrelation of IRAC galaxies for the 
sample described in 13.21 The observed u}{d) is shown in Fig. [5] 
along with the corresponding power-law fit and scaled correlation 
function for DM, calculated as discussed in Fit parameters 
are given in Table [T] The power-law model fits well on the chosen 
scales of 0.3'-10'. The best-fit power law parameters are ro,GG ~ 
3.3 ± 0.3 and 7 = 1.8 ± 0.2, and the best-fit scaled DM model 
yields b% = 2.99 ± 0.40 or bp = 1.73 ± 0.12. 

This accurate value for bp yields bs ~ 3.37 ± 0.82 for the 
S MGs. Converting this t o DM halo mass using the prescription 
of ISheth. Mo & TormerJ ( 1200 ih as described in i]3.4| we arrive at 
log (Mhaio[/i~^ -^o]) = 12.8lQ g. The corresponding halo mass 
for the galaxies is log (Mhaio[/i"^ Mq]) = 11.5 ± 0.2. 

For comparison with other studies that attempted to directly 
measure the autocorrelation function of SMG, it is useful to present 
the SMG clustering in terms of effective power-law parameters for 
their autocorrelation. Assuming linear bias, the SMG autocorrela- 
tion c an be inferred from the cross-correlation by ^ss = ^sg/^gg 
(e.g., ICoiletaL .2009.) . Adopting a fixed 7 = 1.8 for the SMG- 
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Table 1. Correlation results 



Subset 


AT ° 


{z)" 


Power law fif^ 
ro ih~^ Mpc) 7 




Bias model fif* 
bs^G (fe^) {bo) 




Halo mass'^ 

(log/l-l Mq) 


SMGs 
galaxies 


50 
11,241 


2.02 
2.13 


7 7+1.8 
'■'-2.3 

3.3 ±0.3 


1.8 ±0.2 
1.8 ±0.2 


0.8 
1.8 


5.83 ± 1.36 
2.99 ±0.40 


3.37 ±0.82 
1.73 ±0.12 


0.7 
1.8 


12.81^5 
11.5 ± 0.2 



" Number of objects in the SMG sample and in the galaxy sample used for the galaxy autocoiTelation. 
' Median redshift for the SMG sample and for the galaxy sample used for the galaxy autocorrelation. 

Power law model parameters are for the autocorrelation of SMGs (derived from SMG-galaxy projected spatial cross- 
correlation, along with the galaxy angular autocorrelation) and galaxies (derived from their angular autocorrelation). 

Parameters derived from the observed linear fit of the DM model to the observed correlation function, in order to obtain the 
the absolute bias for the SMGs and galaxies (denoted bg and feg, respectively). The linear scaling from the fit corresponds to 
bgbQ for the SMG-galaxy cross-coirelation, and b^ for the galaxy autocon'elation, which in turn yield bQ and bg. 

DM halo mass derived from the absolute bias, using the method described in ii l3.4l 



galaxy cross-correlation, we thus obtain ro,ss ~ ■'^^2% ^ ^ Mpc 
for the autocorrelation of the SMGs. 



4.1 Effects of SMG photo-z errors 

One uncertainty in our estimate of Wp{R) for the SMG-galaxy 
cross-correlation is due to the lack of accurate (that is, spectro- 
scopic) redshifts for roughly half of the SMG population. As de- 
scribed in §[3] in calculating Wp{R) for the cross-correlation, we 
simply assume that the SMGs lie exactly at the best redshifts from 
the photo- 2; analysis of Wardlow et al. (2011). Any uncertainties 
in the SMGs photo-zs could therefore affect the resulting cluster- 
ing measurement. (Note that photo-z uncertainties in the galaxies 
are accounted for implicitly in the correlation analysis, as we uti- 
lize the full galaxy photo-z PDFs.) To examine the effects of SMG 
photo-2 errors, we follow the procedure outlined in § 6.3 of |H1 1| . 
We take advantage of the 44% of SMGs that do have spectroscopic 
redshifts, and determine how errors in those redshifts affect the ob- 
served correlation amplitude. 

Specifically, we shift the redshifts of the spectroscopic SMGs 
by offsets Az/{l+z) selected from a Gaussian random distribution 
with dispersion az/{l + z). To ensure that this step does not arti- 
ficially smear out the redshift distribution beyond the range probed 
by the galaxies, we require that the random redshifts lie between 
1 < z < 3; any random redshift that lies outside this range is dis- 
carded and a new redshift is selected from the random distribution. 
Using these new redshifts we recalculate Wp{R), using the full for- 
malism described in §[3] We perform the calculation 10 times for 
each of several values of 0-^/(1 + z) from 0.05 up to 0.3 (corre- 
sponding to the range of photo-z uncertainties). For each trial we 
obtain the relative bias by calculating the mean ratio of Wp{R), on 
scales 1-10 Mpc, relative to the Wp{R) for the best estimates 
of redshift. We then average the ten trials at each az , and find that at 
most the photo-2 errors cause the clustering amplitude to decrease 
by ~ 10%. The precise magnitude of this effect is unclear given 
the range of uncertainties in the SMG photo-z estimates, but it is is 
significantly smaller than the statistical uncertainties. We therefore 
neglect this effect in our final error estimates. 



4.2 Comparison with previous results 

Here we compare our results to other measurements of SMG clus- 
tering in the literature. The observed clustering ma y depend on the 
flux li mit of the submm sample, as discussed by IWilliams et alj 
JioT l|); measurements of ro that use SMG samples with similar 



submm flux limits are shown in Figure|6ji. Our measurement is sig- 
nificantly more accurate than previous measurements, owing to the 
inclusion of redshift information and the improved statistics in the 
cr oss-correlation. Th e uncertainties are comparable to those quoted 
bv lBlain et all ( |2004|) who estimated ro using counts of close pairs 
in redshift space from spectroscopic surveys. However, these au- 
thors did not account for signifi cant additional sources of error, as 
discussed by lAdelbergeil ( l2005h . Uncertainties in the redshift se- 
lection function for spectroscopic objects, along with the presence 
of redshift spikes and angular clustering of sources, can strongly 
impact the number of expected pair counts for an unclustered dis- 
tribution, and ther efore significantly affect the results for the clus- 
tering ampl itude 1 Adelberget |2005h . In Figure [6ji the large error 
bars for the Isiain et alj 1 2004h point represent the increase in the 
uncertainty by 60% due to a ngular clustering o f sources and red- 
shift spikes (as estimated by lAdelbergej|2005h , but does not in- 
clude the additional uncertainty on the redshift selection function. 
Nonetheless, our measurement of ro is consiste nt with mos t previ- 
ous angular clustering estimates as well as the lBlainet an ( l2004h 
result, and represents a significant improvement in precision. 

As discussed in § 13.41 we convert the observed clustering 
amplitude to A/haio by assuming that SMGs obey simple linear 
bias relative to the dark matter and reside in haloes of similar 
mass. Motivated by the presence of a large overdensity of SMGs 
and powerful s t ar-for ming galaxies in one redshift survey field, 
IChapman etal] ( |2009|) proposed that SMGs obey "complex bias" 
that depends on large-scale environment and merger history, and 
that they may reside in somewhat smaller haloes than would be in- 
ferred from a linear bias model. Future studies using significantly 
larger SMG samples may be able to confirm the existence of more 
complex clustering, but for the present analysis we adopt the sim- 
plest scenario and derive i\fhaio assuming linear bias. 

The characteristic halo mass we measure for SMGs is simi- 
lar to that measured for bright far-lR sources (with fluxes > 30 
mjy at 250 ^m) detected by the Herschel Space Observatory us- 
ing an angular clustering analysis jCoorav e t al. 2010). While it re- 
mains uncertain to what extent bright 250 ^m sources and 850 /^m- 
selected SMGs represent a common population, both samples com- 
prise the luminous end of the star-forming galaxy population de- 
tected at those wavelengths and so may represent physically similar 
systems. In contrast, our obser ved SMG clust e ring i s significantly 
stronger than that reported by lAmblardetai] ( l201lh for "submil- 
limetre galaxies" based on a power-spectrum analysis of Herschel 
350 /im maps, which yields a minimum Mhaio of ~ 3 x 10^^ Mq . 
The differences in clustering amplitude compared to SMGs result 
from the fact that the power spectrum analysis includes unresolved 
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faint sources corresponding to far fainter far-IR luminosities, char- 
acteristic of typical 2 ~ 2 star-forming galaxies rather than the 
powerful, luminous starbursts that are conventionally referred to as 
SMGs in the literature. 



4.3 Progenitors and descendants of SMGs 

Our improved clustering measurement allows us to place SMGs 
in the context of the cosmological history of star formation and 
growth of DM structures. Because the clustering amplitude of dark 
matter haloes and their evolution with redshift are directly predicted 
by simulations and analytic theory, we can use the observed clus- 
tering to connect the SMG populations to their descendants and 
progenitors, estimate lifetimes, and constrain starburst triggering 
mechanisms. 

We first compare the clustering amplitude of SMGs with other 
galaxy populations over a range of redshiftfl Figurej^J) shows the 
approximate ranges of measurements of ro for a variety of galaxy 
and AGN populations. We also show the evolution of ro with red- 
shift for DM haloes of different masses, determined by fitting a 
power law with 7 = 1.8 to the DM correlation function output by 
HALOFIT. Finally, we show the observed ro for the current SMG 
sample, along with the expected evolution in ro for haloes that 
have the observed Mhaio for SMGs at 2: = 2, calculated using 
the median growth rate of haloes as a function of A/haio and z 
jFakhouri, Ma & Bovlan-Kolchinll2010lfl 

Figure |6}? shows that while the DM halo mass for the SMGs 
will increase with time from 2 ~ 2 to 2 = 0, the observed ro 
stays essentially constant, meaning that the progenitors and descen- 
dants of SMGs will be populations with similar clustering ampli- 
tudes. Our measurement of ro shows that the clustering of S MGs 
is co nsistent with optic ally-selected QSOs ( e.g.. Groom et alj 
booi iMvers et al.ll2006l : Ida Angela et alll2008l ; Iross et alj|20m 
SMGs are more strongly clustered than the typical star-formin; 
galaxy populations at all redshifts (e.g. lAdelberger et alj 1200: 
iGilli et al.ll2007l : Iffickox et al.] l2009l : IZehavi et alj|201lh . and are 
clustere d similarly or weaker than massive, passive systems (e.g. , 
Ouadri et ai] l2007l l2008l: l\yake et al.l l2008bl ; iBlanc et al.l l2008l : 
Kim et al.l201 iMZehavi et alj201 ih . The clustering results indicate 



that SMGs will likely evolve into the most massive, luminous early 
type galaxies at low redshift. We note that the descendants of typi- 
cal SMGs are not likely to reside in massive clusters at 2 = 0, but 
into moderate- to high-mass groups of ~a few xlO" h'^ Mq. 
Although some SMGs could evolve into massive cluster galaxies, 
the observed clustering suggests that most will end up in less mas- 
sive systems. 

A schematic picture of the evolution of SMGs is 
shown in Figure [7] which shows evolution in the mass of 
haloes with redshift as traced by their median growth rate 
jFakhouri. Ma & Bovlan-K olchin 2010). The typical progenitors 
of SMGs would have A/haio ~ 10^^ Mq at 2 ~ 5, which cor- 
responds to the host haloes of bright LBGs at those redshifts (e.g., 
iHamana et al] |2004 iLee et al.l l2006h . At low redshift, the SMG 
descendants will have Mhaio = (0.6-5) x 10^^ Mq. Halo 
occupation distribution fits to galaxy clustering sugg est that these 
haloes host galaxies with luminosities L ~ 2-31/* jZehavi et alj 



1 iMvers et al] )20Qd) and lRoss et all )2009h determine ro from QSOs as- 
suming a power law correlation function with 7 = 2. To estimate ro for 
7 = 1.8, we multiply the quoted values by 0.8, appropriate for fits over the 
range 1 < i?. < 100 h'^ Mpc. 
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Figure 6. (a) Our new measurement of the autocorrelation length ro for 
SMGs, compared to previous results using sam ples with similar ~850 /^m 
flux limits. The two sets of error bars on the iWebb et all )2003h measure- 
ment indicate statistical (±3 Mpc) and sys tematic (±3 Mpc) un- 
certainties separately. On the iBlain et all j2004l) measurement, the smaller 
eiTors represent the uncertainties quoted by the authors, while the larger 
errors account for angular clustering and redshift spikes as estimated by 
lAdelberge 3 i2005l) . Our results are consistent with previous measurements 
and represent a significant improvement in precision, (b) Our measure- 
ment of the autocorrelation length ro of SMGs, compared to the approxi- 
mate ro (with associated measurement uncertainties) for a variety of galaxy 
and AGN p opulations: optically-selected SDSS QSOs at < 2 < 3 
( IMvers et all 2006; Ross et al. 2009), Lyman-break galaxies (LBGs) at 
1.5 < 2 < 3.5 (.Adelberger et_ al. 2005), MIPS 24 |j.m-selected star-forming 
galaxies at < z < 1.4 loilli et al. 20(3), typical red an d blue galax- 
ies at 0.25 < z < 1 from the AGES Mickox et al.ll2009h and DEEP2 
et al.ll200^ spectroscopic s urveys, luminous red galaxies (LRGs) at 
< 2 < 0.7 I'wake et al' 2008b), and optically-selected galaxy clusters at 
0.1 < 2 < 0.3 (Estrada, Sefusatti & Frieman 2009). In addition, we show 
the full range of ro for low-redshift galaxies with r-band luminosities in the 
range 1.5 to 3.5 L*, derived from the luminosity dependence of clustering 
presented by Zehavi et al. (201 1); these luminous galaxies are primarily el- 
lipticals, as discussed in ^ 14.31 Dotted lines show ro versus redshift for DM 
haloes of different masses. The thick solid line shows the expected evolu- 
tion in ro , accounting for the increase in mass of the halo, for a halo with 
mass corresponding to the best-fit estimate for SMGs at 2 = 2. The results 
indicate that SMGs are clustered similarly to QSOs at 2 ~ 2 and can be 
expected to evolve into luminous elliptical galaxies in the local Universe. 
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Figure 7. Broad schematic for the evolution of halo mass versus redshift 
for SMGs, showing the approximate halo masses corresponding to likely 
progenitors and descendants of SMGs. Lines indicate the median growth 
rates of haloes with redshift jpakhouri. Ma & Bovlan-KolchiijlioTol) . SMG 
host haloes are similar to those those of QSOs at 2 2, and correspond to 
bright LBGs at ^ ~ 5 I'Hamana et al."2004";'Lee et al."2006^ and ~ 2-3L* 
ellipticals at z = jzehavi et al. 201 1; Stott et al. 201 1) . 



a population dominate d by ellipticals with predominantly 
slow- rotating kinematics (e.g. jTempel et alj201 iLICappellari et alJ 
1201 ih . Assuming typical mass-to-lig h t ratio s for massive galaxies 
(e.g., ISaldry, Glazebrook & Driveii 1200 Sl) , these luminosities 
correspond to stellar masses -(1.5-2.5) X 1O^^M0, in close 
agreement with direct measurements of the relationship between 
halo mass and central galaxy stellar mass for X-ray selected groups 
and c lusters, for which log M* x 0.27 log A/haio + 7.6 dStott et alj 
l201lh . 



4.4 SMG lifetime and star formation history 

We next estimate the SMG lifetime, making the simple assumption 
that every dark matter halo of similar mass passes through an SMG 
phas^E so that 



tSMG = At 



n-SMG 
Mhalo ' 



(8) 



where At is the time interval over the redshift range covered by the 
SMG sample, and tismg and Whaio are the space densities of SMGs 

and DM haloes, respectively. 

Using the halo mass function of [Tinker et"al ] tOOA the 
space density of haloes with log (Mhaio[ft~^ Mq]) = 12.8to,l 
is dnhaio/dlnM = {2.1^11) x 10"'' Mpc"^. We adopt a space 
density of SMGs at z ~ 2 of ~ 2 x 10~ ^ Mpc~^, correspond- 
ing to results from previous surveys (e.g.. IChapman et"^l2005l : 
ICopDinetalJl2006l : Schael et al. in preparation). This density is 
~ 50% higher than that observed in the LESS field IWardlow et alj 



^ Note that here we use the median growth rate of haloes, which for haloes 
of ~ 10^3 h-^ Mq is Ri35% lower than the mean growth rate, owing to 
the long high-mass tail in the halo mass distribution. 

^ If the average halo experiences more or fewer SMG phases in the given 
time interval, the lifetime of each episode will be correspondingly shorter 
or longer, respectively. 



l201lh . which has been shown to contain a somewhat smaller den- 
sity of SMGs compared to other surveys (WeiB et al. '2009'). 

The ratio of these space densities yields a duty cycle (the frac- 
tion of haloes that host an SMG at any given time) of ~ 10%. We 
assume the SMGs occupy the redshift range 1.5 < z < 2.5 , which 
includes roughly half of the SMGs in the IWardlow etal] j201lh 
sample and corresponds to At = 1.6 Gyr. We thus obtain a life- 
time for SMGs of tsMG = 110tgo° Myr. Clearly, even our im- 
proved measurement of SMG clustering yields only a weak con- 
straint on the lifetime, but this is consistent with lifetimes esti- 
mated ^xjm gas consump tion times and star-formation timescales 
(e.g., lGreveetalJl2005l ; iTacconi et al.ll2006l ; iHainline et all 1201 ih 
and theoretical models of SMG fueling through mergers (e.g. , 
Mihos & Hemauistlll994l : ISpringel, Pi Matteo & Hemquisil2005l : 



Narayanan et alj|2010l) 



Constraints on SMG descendants from clustering can also 
yield insights into their their formation histories. Measurements 
of the stellar plus molecular gas masses of SMGs from SED fit- 
ting and dynamical studies are in the range ~ (1-5) X 10" Mf^ 
JSwinbank et al. 2006; Wardlow et al. 201lJ; IHainline et a"i]|201l| - 
llyison et al. 201 1 ; Michalowski et al. 201 IJ. While these estimates 
can be uncertain by factors of a few, they are in a similar range to 
the stellar masses of SMG descendants as indicated by their clus- 
tering, as discussed above. This correspondence suggests that if a 
significant fraction of the molecular gas is converted to stars dur- 
ing the SMG phase, then these galaxies will subsequently experi- 
ence relatively little growth in mass from z ~ 2 to the present. 
This in turn puts limits on the star formation history. Star-forming 
galaxies at jz ~ 2 ty pically exhibit spe cific star formation rates of 
M*/M* ~ 2 Gyr-i l lElbaz et alj2"oill) . at which the SMGs would 
only need to form stars for 500 Myr in order to double in mass. 
We may therefore conclude, from the clustering and stellar masses 
alone, that the SMGs evolve from star-forming to passive states 
relatively quickly (within a Gyr or so) after the starburst phase, 
and that the descendants spend most of their remaining time as 
relatively passive systems. This scenario is consistent with mea- 
surements of the stellar populations in ~ 2-3 L* ellipticals, which 
have typical ages of ~10 G yr and show little evidence for younger 
components (e.g., |Nelan e t al. 2005; Allanson etal. 2009), imply- 
ing that the vast majority of stars were formed above z ~ 2 with 
little additional star formation at lower redshifts. 

The halo masses of SMGs may also provide insight into 
the processes that prevent their descendants from forming new 
stars. Star formation can be shut off rapidly at the end of the 
SMG phase, either by ex haustion of the gas supply, or by energy 
input from a QSO (e.g., [P i Matteo, Springel & HernquistI l200a ; 
ISpringel, Pi Matteo & Hern quist 2005). Powe rful winds are ob- 
served in liuninous AGN (e.g., f eruglio e t alj|2010l ; iFischer et alj 
l20ld ; ISturmetal.ll201ll ; iGreene et al4,20llh an d have also been 
seen m some SMGs (e.g.. lAlexander et alj|20ld Harrison et al. in 
preparation), although for the SMGs is unclear whether the winds 
are driven by the starburst or AGN. Even if the formation of stars 
is rapidly quenched, over longer timescales the galaxy would be 
expected to accrete further gas from the surro unding halo, result- 
ing in significant a dditional star formation (e.g., lBower et alj2006l ; 
ICrotonetal]|200 6). Recent work suggests that energy from ac- 
creting supermassive black holes, primarily in the form of radio- 
bright relativistic jets, can couple to the hot gas in the surround- 
ing halo, producing a feedback cycle that prevents rapid cooling 
(e.g., lRaffertv, McNamara & Nulsenll2008l) . This mechanical black 
hole feedback is an key ingr edient of s uccessful models for the pas- 
sive galaxy population (e.g., ICroton et aL,2006, ; ,Bower et al.i200^ ; 
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Figure 8. Redshift distribution of LESS SMGs I'War dlow et al.l201 ih . coin- 
pared to the simple models for SMG triggering based on the rate at which 
haloes cross a threshold mass Mthrosh = 6 X 10^^ Mq (see §|4!5). 
The uncertainties i n the number counts are an approximation of Poisson 
counting statistics <Gehrels|[T983) . The black dotted line shows the (arbi- 
trarily normalized) number of haloes crossing this threshold in each red- 
shift interval (Equation |9) while the dashed red line shows this distribution 
multiplied by the evolution in the molecular gas fraction (Equation IIOK 
where /mol is taken from the model predictions of Lagos et al. (2011) and 
is shown by the gray dot-dashed line. The remarkable agreement between 
the second model and the observed number counts suggests that the evolu- 
tion of the SMG population can be described simply in terms of two quan- 
tities: the growth of DM structures and the variation with redshift of the 
molecular gas fraction in galaxies. 



iBower. McCarthy & Bensoij|2008l; ISomerville et al.ll200^ . Inter- 
estingly, the clustering of radio galaxies at 2 < 0.8 in dicates that 
they re side in haloes of mass > 10" h-^ Mr^ (e g.. IWake et al 
20Q8a: 'Hickox et al."2009l: iMandelbaum et alj|2009l : lDonoso et" 
2010; Fine et al. 2011), precisely the enyironments that will host 
the descendants of SMGs. Thus the strong observed clustering for 
SMGs can relate them directly to the radio-bright active galactic 
nucleus population that may regulate their subsequent star forma- 
tion. 



4.5 Evolutionary links with QSOs and the SMG redshift 
distribution 

Finally, the observed clustering of SMGs provides insights into 
the processes that trigger and (possibly) shut off their rapid star 
formation activity. As discussed in § [T] powerful local starbursts 
(i.e. ULIRGs) are predominantly associated with major mergers 
and appear to be associated with the f ueling of luminous Q SOs as 
part of an evolutionary sequence (e.g.. lSanders et af]| 19881) . How- 
ever it is unclear if a similar connection exists between SMGs and 
high-2 QSOs. One robust prediction of any evolutionary picture is 
that SMGs and QSOs must display comparable large-scale clus- 
tering, since the evolutionary timescales are significantly smaller 
than those for the growth of DM haloes. At all redshifts, QSOs 
are found in haloes of similar mass ~ a few xlO^^ h-^Ma^ 



are lounu m naioes or similar mass ~ a lew x iu n iVjQ 
(e.g..lCroometalJl2005l ; iMvers et aUbood Ida Angela et"ai]|2008l : 
IRoss et alj l2009l : Figure |6). The characteristic Mhaio provides 
a strong constraint on models of Q SO fueling by the major 
mergers of gas-rich galaxies (e.g.. iKauffmann & Haehneltll200(]| : 



ular instabilities (e.g.. lMo. Mao & Whitelll998l;lBower et albood 
lOenze l et al. 20 ^ or accreti o n of recycled cold gas from ev olved 
stars jCiotti & Ostrikeill2007l : ICiotti. Ostriker & ProgjbOlOl) . and 
is similar to the mass at whic h galaxy populat i ons transition from 
star-forming to passive (e . g.. ICoil et al.ll2008l: iBrown et all |2008|; 
IConrov & Wechslej |2009| ; iTinker & Wetzelll2010l) . The observed 



clustering of SMGs at 2 ~ 2 from the present work is con- 
sistent with that for QSOs, as well as hig hly active obscured 
objec ts including powerful obscure d AGN |H11|; AUevato et al] 



l201lh and dust-obscured galaxies jBrodwin et al 



2008h . Thus 



these may indeed represent different phases in the same evolu- 
tionary sequence, and energy input from the QSO may be re- 
sponsible for the rapid q uenching of star formation at the end 
of the SMG phase (e.g., [P i Matteo, Springel & Hernauistll2005l ; 
ISpringel. Pi Matteo & Her nquist 200 as discussed in ^14.41 

A connection with QSOs may imply that triggering of SMGs 
is also related (at least indirectly) to the mass of the parent PM 
halo. In this case, the evolution of large-scale structure may broadly 
explain why the SMG population peaks at 2 ~ 2.5 and falls at 
higher and lower redshifts. In the simplest possible such scenario, 
SMG activity is triggered when the halo reaches a certain mass 
-A/halo ~ Mthrosh (sec Figure 16 of iHickox et alj l2009i for a 
schematic illustration of this picture). In a given volume, the num- 
ber of haloes crossing this mass threshold as a function of redshift 
is: 

!±££!!i (X nhalo(Afthrcsh, 2)Mhalo(Mthrcsh, z)tsMG—;—, (9) 

dz az 



where Whaio and Mhaio are the n umber density (e.g., lTinker ( 



2008 ) and typical growth rate jFakhouri. Ma & Bov lan-Kolchin 
2010l) . respectively, of haloes of mass Mthresh at redshift z, tsMG 



is the SMG lifetime, and dV/dz is the differential comoving vol- 
ume over the survey area. If an SMG is triggered every time a halo 
reaches A/throsh, then the observed number density of SMGs will 
be proportional to dA^thrcsh/dz. However, the huge star forma- 
tion rates of SM Gs require a larg e reservoir of molecular gas (e.g., 
fGreve et al, 200 ji lTacconi et alj2006 . 2008). a nd the molecular gas 
fraction increases strongly with reds hift (e.g., lTacconi et al]|201(]| : 
iGeach et alJlioTTI : iLagos et alJlioTlh . This evolution may explain 
why the most powerful starbursts at low re dshift (ULIRGs) have 
lower typical SFRs tha n z ~ 2 SMGs (e.g.. lLe Floc'h et al.|[20o3 : 
lRodighieroetalJl2010l) . Therefore it may be reasonable to assume 
that the number counts of SMGs also depend on fmoi, with the 
simplest possible prescription being: 

J « ^ /mol 2 . (10) 

dz dz 
In Figure [8] we show the ob served redshift distribution of 
LESS SMGs (W ardlow et alj|201lh . compared to the distributions 
predicted by Equations ^ and ilOi . assuming Mthresh = 6 x 10^^ 
h^^ AIq. For simplicity, the evolution in /moi i s taken from pre- 
dictions of the GALFORM model of iLagos et 'al. (2011), which 
agree s broadly with observations (see Figure 2 of IGeach et alj 
I2OI ll) and so provides a simple parameterisation of the current em- 
pirical limits on the molecular gas fraction in galaxies. It is clear 
from Figure[8]that there is remarkable correspondence between our 
extremely simple prescription and the observed redshifts of SMGs. 
Of course this "model" does not account for a wide range of pos- 
sible complications and the normalisations of the distributions are 
arbitrary. However, this exercise clearly demonstrates that if SMGs, 
like QSOs, are found in haloes of a characteristic mass, then their 
observed redshift distribution may be explained simply by two ef- 
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fects: the cosmological growth of structure combined with the evo- 
lution of the molecular gas fraction. Thus SMGs likely represent a 
short-lived but universal phase in massive galaxy evolution, associ- 
ated with the transition between cold gas-rich, star-forming galax- 
ies and passively evolving systems. 



5 CONCLUSIONS 

In this paper we measure the cross-correlation between SMGs and 
galaxies in the LESS survey of ECDFS, and observe significant 
clustering at the > 4a level. We obtain an autocorrelation length for 
the SMGs of ro = 7-7^2,1 h'^ Mpc, assuming 7 = 1.8. This clus- 
tering amplitude corresponds to a characteristic DM halo mass of 
log (Mhaio[/i"^ Mq\) = 12.810 5. Using this estimate of A/haio 
and the space density of SMGs, we obtain a typical SMG lifetime 
ofisMG = 11018^° Myr. 

The observed clustering indicates that the low-redshift descen- 
dants of typical SMGs are massive (~2-3 L*) elliptical galax- 
ies at the centers of moderate- to high-mass groups. This predic- 
tion is consistent with pre vious suggestions based on the dynami- 
cal dSwinbank eTai]|2006l) and stellar masses (e.g.. ^Hainline et aT] 
l201lh of SMGs, and is also consistent with observations of local 
massive ellipticals, which indicate that they formed the bulk of their 
stars at z > 2 and have been largely passive since. The clustering 
of SMGs is very similar to that observed for QSOs at the same 
redshifts, consistent with evolutionary scenarios in which SMGs 
and QSOs are triggered by a common mechanism. Assuming that 
SMGs, like QSOs, are transient phenomena that are observed in 
haloes of similar mass at all redshifts, the redshift distribution of 
SMGs can be explained remarkably well by the combination of the 
cosmological growth of structure and the evolution of the molecu- 
lar gas fraction in galaxies. 

This accurate clustering measurement thus provides a valuable 
observational constraint on the role of SMGs in the cosmic evolu- 
tion of galaxies and large-scale structures. We conclude that SMGs 
likely represent a short-lived but universal phase in massive galaxy 
evolution that is associated with the rapid growth of black holes 
as luminous QSOs, and corresponds to the transition between cold 
gas-rich, star-forming galaxies and passively evolving systems. 
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